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Abstract 

The production of the X(3872) is studied in pp colHsions at -\/s = 7TeV, using decays 
to j/xpn^n^ , where the ]/xp decays to two muons. The data were recorded by the 
CMS experiment and correspond to an integrated luminosity of 4.8fb^\ The mea- 
surements are performed in a kinematic range in which the X(3872) candidates have 
a transverse momentum 10 < px < 50GeV and rapidity \y\ < 1.2. The ratio of 
the X(3872) and tp{2S) cross sections times their branching fractions into j/xpn^n^ is 
measured as a function of pj. In addition, the fraction of X(3872) originating from 
B decays is determined. From these measurements the prompt X(3872) differential 
cross section times branching fraction as a function of pj is extracted. The k^tt^ 
mass spectrum of the ^/xpn^n^ system in the X(3872) decays is also investigated. 
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1 Introduction 

The discovery of the X(3872) resonance by the Belle experiment in 2003 111 and its subsequent 
confirmation by BaBar, CDF, and DO lUfll has attracted a large interest in "exotic" quarko- 
nium spectroscopy since it was the first observation of an unexpected charmonium candidate. 
Many new unconventional states with masses above the open-charm or open-bottom thresh- 
olds, m > m(DD)andm > m(BB),respectively, have been observed [5 J. There are several inter- 
pretations of the X(3872) state: a charmonium state, a D*D molecule, or a tetra quark state |5|. 
The X(3872) has been observed in several decay channels, including ]/tpn^n^ , D*D, J/t/'7, 
i/?(2S)7, and ]/tpLO. The analysis of the X(3872) angular distributions in decays to ^/xpn^n^ 
favours J^'~ = 1++ or 2 ^ IISIIZI- The inclusive production cross section of the X(3872) res- 
onance has been measured by the LHCb experiment |8l. At the Tevatron it was observed 
that the X(3872) is produced both through "prompt" processes, in which the X(3872) reso- 
nance is created directly, and through decays of B hadrons |9|, generally referred to as "non- 
prompt". Experimentally, nonprompt processes are distinguishable through the displacement 
of the X(3872) decay vertex from the primary vertex. Prompt production of quarkonium states 
in proton-proton collisions is usually described in the framework of nonrelativistic quantum 
chromodynamics (NRQCD) [5J. Quantitative predictions have been calculated for the differ- 
ential production cross section of the X(3872) in pp collisions at the Tevatron and pp collisions 
at the Large Hadron Collider (LHC) |TOl [TT| . Measurement of the prompt production rate at 
the LHC as a function of transverse momentum provides a test of the NRQCD factorization 
approach to X(3872) production. 

In this paper a measurement of the differential X(3872) production cross section is presented 
using decays into ^/xpn^n^, with the subsequent decay of the ]/xp into a pair of muons. The 
analysis makes use of pp collision data recorded by the Compact Muon Solenoid (CMS) exper- 
iment at the LHC in 2011, at a centre-of-mass energy of 7TeV, corresponding to an integrated 
luminosity of 4.8fb^\ The analysis is performed in the kinematic range of pj of the j/xpn^TZ^ 
system between 10 and 50GeV and the rapidity within \y\ < 1.2. The cross section measure- 
ment proceeds by determining the ratio of the X(3872) and J/'(2S) cross sections, where both 
states decay to ]/ipTi^7z^. In this ratio, systematic uncertainties common to both states largely 
cancel, either partially, as those related to the trigger and the reconstruction of the ]/tp mesons, 
or fully, as for the integrated luminosity. The fraction of nonpromptly produced X(3872) states 
is also measured. The cross section times branching fraction for prompt X(3872) production 
with ]/xpTT'^n~ in the final state is then extracted by using a previous CMS measurement of the 
differential cross section for prompt ^{2S) production in the same kinematic range [12J. The 
differential cross section for prompt X(3872) production times the branching fraction is deter- 
mined for the first time as a function of transverse momentum. Finally, the invariant-mass 
distribution of the dipion system in X(3872) ]/ipn^n^ decays is studied. 

Throughout the analysis, the corrections for detector acceptances and efficiencies are deter- 
mined under the assumption that the X(3872) has quantum numbers J^*-- = 1++ and that both 
the X(3872) and the ^{2S) are unpolarized. The unknown polarizations of the X(3872) and the 
j/7(2S) lead to large uncertainties, in particular in the acceptance of the final-state muon pair for 
extreme polarization hypotheses. 

This paper is structured as follows: after a brief description of the CMS detector in Section |2| 
the data sample and event selection are discussed in Section |3j The measurement of the cross 
section ratio is reported in Section |4j Section [5| gives the measurement of the relative fraction of 
nonprompt X(3872) production. In Section |6|the cross section for prompt X(3872) production 
is presented. Finally, in Section ItI the tt+tt^ mass spectrum in X(3872) to ]/ipn^W decays is 
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reported. 

2 CMS detector 

The central feature of the CMS apparatus is a superconducting solenoid of 6 m internal di- 
ameter. Within the solenoid, in a 3.8 T magnetic field, are a silicon pixel and strip tracker, a 
lead tungstate crystal electromagnetic calorimeter, and a brass /scintillator hadron calorime- 
ter. The main subdetectors used in this analysis are the silicon tracker and the muon sys- 
tem. Charged-particle trajectories are measured in the pseudorapidity region \tj\ < 2.5, where 
y] = — ln[tan(0/2)], with 9 being the polar angle with respect to the anticlockwise-beam direc- 
tion. The tracker provides an impact parameter resolution of ?»15/^m. Muons are detected in 
three types of gas-ionization detectors embedded in the steel flux return yoke: drift tubes in the 
barrel, cathode strip chambers in the endcaps, and resistive plate chambers in both the barrel 
and endcaps. Matching the muons to the tracks measured in the silicon tracker results in a 
transverse momentum resolution between 1 and 1.5%, for px values up to 50 GeV. A two-level 
trigger system selects relevant pp collision events for offline reconstruction. The first level (LI) 
of the CMS trigger system is composed of custom hardware processors. The LI trigger condi- 
tions are adjusted such as to limit the trigger rate to less than 100 kHz. The high-level trigger 
(HLT) runs on a processor farm to further reduce the rate to a few 100 Hz before data storage. 
A detailed description of the detector can be found elsewhere llT3l . 

3 Event selection 

The event selection criteria are largely driven by requirements imposed at the trigger level. At 
both trigger levels, at least two muons are required. At the HLT, events are accepted if the two 
muons are of opposite charge, have an invariant mass between 2.95 and 3.25 GeV, vertex fit 
probability greater than 0.5%, and rapidity < 1.25. In 2011, the transverse mo- 

mentum threshold at the trigger level for the dimuon system was initially 6.9 GeV, which was 
increased to 9.9 GeV near the end of data taking. In addition, to cope with increasing instanta- 
neous luminosities, events in which two muons bend toward each other in the magnetic field 
were rejected by criteria added near the beginning of the data taking. The data sample consists 
of events where an average of six pp collisions in the same bunch crossing (pileup) occur. 

In the offline event selection, similar criteria are imposed on the muon pair. Muons are required 
to have opposite sign. The rapidity of the muon pair is required to be \y{}i^}i^)\ < 1.25. 
For the first part of the data, 2011a, corresponding to an integrated luminosity of 2.1 fb^^, a 
minimum dimuon transverse momenta of pT{fi^ji^) = 7 GeV is required. For the second 
part of the data, 2011b, corresponding to an integrated luminosity of 2.7 fb^^, the transverse 
momentum threshold is increased to Pt(h^M ) — 10 GeV. The muon identification criteria 
are very similar to those used in a previous CMS analysis [12J. Each candidate muon track 
must be matched to a triggered muon and have a transverse momentum pT(f<) > 4 GeV in 
the central-pseudorapidity interval \y|{^)\ < 1.2, or Pt{}i) > 3.3 GeV in the forward region 
1.2 < < 2.4. These requirements, together with the rapidity and transverse momentum 

selection criteria on the muon pair, define the ]/xp acceptance, AiJ/ip). Each muon track must 
have at least 11 tracker hits, of which at least two are in the silicon pixel layers. The tracks 
are required to intersect the beam line within a cylinder of 3 cm in radius and 30 cm in length 
around the primary vertex position, selected as the vertex with the largest sum of pj of the 
tracks associated with it. The track fit is required to have a per degree of freedom, x^/ ndf, 
smaller than 1.8. The dimuon vertex fit probability is required to be above 1%. The invariant 
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mass of the muon pair is required to be within 75 MeV of the fitted ]/xp peak, corresponding to 
about ±2.5 times the detector resolution for the ]/xp mass region. Based on this selection, an 
almost background-free sample of about ten million reconstructed ]/tp candidates is obtained. 

The ]/xpK^n^ system is reconstructed by combining the candidate muon tracks from each 
candidate ]/tp with pairs of oppositely charged tracks, which are assumed to be pions. Each 
ji^ji^n^n^ combination is refitted, constraining the four tracks to come from a common ver- 
tex and the muon-pair invariant mass to the ]/xp mass ||14|. Combinations yielding a vertex 
fit probability smaller than 5% are rejected to suppress combinatorial background. The pion 
tracks must have a fit x^/ ndf < 5 and contain at least two (seven) silicon pixel (strip) hits. The 
refitted pion tracks must also have a transverse momentum larger than 600 MeV. 

Random combinations of tracks form a significant combinatorial background to the X(3872) 
and tp{2S) signals. For signal events, the pions are expected to have a direction close to that 
of the j/ip candidate. Exploiting this property, the combinatorial background is reduced by 
requiring the distance AR = ^ {Ar]Y + {A(pY, where A?/ and are the pseudorapidity and 
azimuthal angle differences between the pion and the ]/'ip candidate momenta, to be smaller 
than 0.55. The requirement AR < 0.55, together with the pion transverse momentum selection, 
define the dipion acceptance, A(7r+7r~). The event selection criteria are driven by studies 
from simulation, whose description is reported below, using the quantity S/i/S + B, where 
S and B represent the numbers of signal and background candidates, respectively, in a ±2(7" 
window around the X(3872) mass, and a is the mass resolution (about 6 MeV). In addition, 
the Q value of the decay is required to be smaller than 300 MeV, where Q = m[}i^ ji^ n^) — 
niQ/tp)^^^ — m( 7r+ TT^ ), with m(f/+f/^ 71+ TT^) being the invariant mass of the J/t/^/r+Tr^ system, 
m[7i^TZ^) the invariant mass of the pion pair, and m(J/i/;)™*^ the world-average ]/tp mass ||T4| . 
This selection criterion constrains the mass of pion pairs from X(3872) decays to values larger 
than about 470 MeV and removes about 20% of the remaining background, while retaining 
97% of the X(3872) signal, as determined from simulation. The invariant mass of the pion pair 
measured in data, as shown in Section]^ has a negligible contribution below 500 MeV, and thus 
no bias is introduced with the Q-value requirement. 

The ]/tpn^n^ candidate is required to be in the rapidity region \y\ < 1.2, to have transverse 
momentum pj < 50GeV and pj > lOGeV for the period 2011a or pj > 13.5 GeV for the 
period 2011b. The resulting data sample consists of about 1.9 million ]/ipTT'^n^ candidates 
with invariant mass between 3.6 and 4.0 GeV. The average number of ]/ipn^n^ candidates 
per event for events with at least one such candidate is reduced from 7.8 to 2.2 after the event 
selection. 

Figure [l] shows the invariant-mass distribution of the ]/ip7T'^7i~ candidates passing the full 
event selection. An unbinned maximum-log-likelihood fit is used, where Gaussian distribu- 
tions describe the signals and a Chebyshev polynomial the background. More details about 
the fit will follow in Section [i] Clear ^{2S) and X(3872) signals are observed with widths of 
about 5 MeV and 6 MeV, respectively, dominated by the detector resolution and consistent with 
simulation. 

Detailed event simulations are used to determine detector effects such as acceptances, efficien- 
cies, and resolutions. Events containing X(3872) or tp{2S) states are generated using pythia IflSlI , 
and decayed using evtgen [16|, with the signal resonances forced to decay into the j/xpn^Ti' 
final state. Photon final-state radiation (FSR) is implemented using PHOTOS ITtI ITsI . The 
X(3872) and tp[2S) resonances are assumed to be unpolarized. Since pythia does not include 
the simulation of X(3872) production and decay, the program is modified to use the Xci particle 
with its mass set to 3871.6 MeV. The Xci particle has the quantum numbers J^*-- = 1++, corre- 
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Figure 1: The J/i/?7r+ tt^ invariant-mass spectrum for 10 < px < 50 GeV and \y\ < 1.2. The 
Hnes represent the signal-plus-background fits (solid), the background-only (dashed), and the 
signal-only (dotted) components. The inset shows an enlargement of the X(3872) mass region. 



sponding to those favoured for the X(3872) |ISl[T9l. Simulated events for prompt production 
are used as the baseline. Events with B-hadron decays are simulated and used in the X(3872) 
nonprompt-fraction measurement. The X(3872) — > j/xpn^n^ decay is generated with an in- 
termediate resonance, as suggested by previous measurements ||Zl|20| and confirmed in this 
analysis (Section|7]|. In evtgen a two-body phase-space decay is used for the X(3872) ]/ipp^ 
decay, and the decay to a pair of pions is generated with decay-angle distributions reflecting 
their respective spins. A nonresonant X(3872) — >• j/tpn^n^ decay is also considered using the 
EVTGEN model for the tp{2S) — ?> j/tpn^n^ decay. The study of systematic uncertainties uses 
a version of PYTHIA that includes colour-octet contributions with NRQCD matrix elements, as 
determined from CDF data 112111221. 

Large samples of simulated events are produced separately for the X(3872) and tp{2S) res- 
onances, both for prompt production and nonprompt production in B-hadron decays. The 
response of the detector is simulated in detail using Geant4 [23|. The simulated samples 
are processed through the trigger emulation and event reconstruction of the CMS experiment, 
without taking into account other pp collisions in the same bunch crossing (pileup) since the 
analysis is not sensitive to it, as discussed in Section |4] 

4 Measurement of the cross section ratio 

The ratio of the cross section times the ^/xpn^n^ branching fraction is obtained from the mea- 
sured numbers of signal events for X(3872) and ip{2S), Nx(^3S72) ^rid N^(2S)/ correcting for the 
efficiency (e) and acceptance (A) estimated from simulations, according to 



tr(pp ^ X(3872) + anything) ■ g(X(3872) ^ J/t/^Ti+Ti-) _ Nx(3872) " ^i/,(2S) " g|/>(2S) 
C7(pp ^ ip{2S) + anything) ■ B{ip{2S) ]/\pn+n-) ~ N^^2S) " ^x(3872) " ex(3872) ' 

(1) 
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The acceptance corrections account for the kinematic reach of the dimuon trigger and the an- 
gular acceptance of the CMS detector. These corrections depend on assumptions about the 
angular distribution of the final-state muon and pion pairs. To minimize the effect of these 
assumptions, the measurement is also presented as a "fiducial" cross section ratio, defined as 



R 



fiducial 



^X(3872) ■ ^i/;(2S) 
^!/)(2S) ■ £x(3872) 



(2) 



within a phase-space window with the following kinematic requirements on the muons, dimuons, 
andpions: muons with pT(f<) > 4GeV for | < l.Zandpxlf^) > 3.3GeV for 1.2 < \f]{j^)\ < 
2.4; pj{ji^fi^) > 7GeV and < 1.25 for the dimuons; each pion with transverse 

momentum greater than 600 MeV and a distance with respect to the dimuon AR < 0.55. 

The signal yields are determined from unbinned maximum-likelihood fits to the invariant-mass 
spectra of the ]/xpn^TZ^ system, separately for the X(3872) and xp{2S), in the mass windows 
3.75-4 GeV and 3.6-3.8 GeV, respectively, and in five bins of pj with edges: 10, 13.5, 15, 18, 
30, and 50 GeV. Following the evolution of the trigger thresholds with time, the first bin in 
transverse momentum, 10-13.5 GeV, includes only data from the period 2011a, while for pj 
bins above 13.5 GeV, the full dataset (2011a+2011b) is used. The inclusive signal yield for pj 
between 10 and 50 GeV is determined by combining the first pj bin from 2011a, weighted to 
account for luminosity and trigger differences, with the remaining bins from the full dataset. 

In the fits, the ^{2S) resonance shape is parametrized using two Gaussian functions with a 
common mean, while a single Gaussian is used for the X(3872) signal. The nonresonant back- 
ground is fitted with a second-order Chebyshev polynomial. The free parameters in the fit are 
the signal and background yields, the mass and widths of the Gaussian functions, the frac- 
tion of signal associated with each Gaussian, and two background-shape parameters. Figure |2] 
shows examples of fitted mass distributions for a low- and a high-transverse-momentum bin. 
The measured numbers of X(3872) and ip{2S) signal events are listed in Tablejl] 

The acceptances and efficiencies of the X(3872) and tp{2S) final states are factorized into four 
components, each of which is determined individually from the simulation: the acceptance 
A{]/tp) and efficiency e(J/j/') for the trigger and detection of the ]/tp, and the acceptance A ( 7r+ tt^ ) 
and efficiency e(7r^7r^) for the pion pair, including the ^/xpTz^n^ vertex probability require- 
ment. The acceptances are the same for the 2011a and 2011b datasets for the pj bins in common 
(px > 13.5 GeV). The efficiency is calculated for the 2011a dataset in each bin since the changes 
in efficiency related to the trigger evolution during data taking do not affect the efficiency ratio. 
The average value of A ■ e in each pj bin is determined using fine-grained bins in transverse 
momentum as 

Table 1: Measured numbers of signal events, Nx{3S72) N^(2S)/ ^rid the ratios of the X(3872) 
and i/'(2S) efficiencies (e) and acceptances (A) as a function of the ]/ipn^7i^ pj. For the first 
transverse momentum bin only the data from period 2011a are included. All uncertainties are 
statistical only. 



Dataset 


Pt (GeV) 


^X(3872) 


^i/>(2S) 






£X(3872) 


^X(3872)''=X(3872) 


2011a 


10-13.5 


1850 ± 200 


25 450 ± 330 


1.055 ±0.011 


0.999 ± 0.025 


2011a+b 


13.5-15 


1700 ± 170 


24 130 ± 440 


1.032 ± 0.014 


0.951 ± 0.025 


2011a+b 


15-18 


2770 ± 210 


39 450 ± 470 


1.031 ± 0.011 


0.979 ± 0.020 


2011a+b 


18-30 


3360 ± 230 


56 920 ± 510 


1.035 ± 0.011 


1.019 ±0.018 


2011a+b 


30-50 


860 ± 140 


12 130 ± 230 


1.052 ± 0.037 


1.103 ± 0.056 


2011a+b 


10-50 


11 910 ±490 


178 540 ± 850 


1.040 ± 0.006 


0.984 ± 0.017 



6 



4 Measurement of the cross section ratio 




1500 



3000. 



3.75 



3.8 



3.85 



3.9 3.95 A 

m(J/\|/ TtV) [GeV] 



4 



3.75 



3.8 



3.85 



3.9 3.95 4 

m(J/v|/ 71 V) [GeV] 



4 



Figure 2: The j/xpn^n^ invariant-mass distribution in the X(3872) region for two bins of trans- 
verse momentum, 10-13.5 GeV (left) and 18-30 GeV (right). The lines represent the signal-plus- 
background fits (solid) and the background-only components (dashed). The ;t^/ndf of the fit 
is also reported. 



where Ni is the number of signal events observed in the data. A' = A'Q/xp) ■ A^n^n^), e' = 
e'{]/tp) ■ e'{TZ^TT^) are the acceptance and efficiency in each fine bin, and N^^^ is the number 
of fine bins contained in each pj interval. This procedure accounts for the large variation in 
acceptance and efficiency over the wide pj bins, relying on the pj spectrum from the data. The 
number of signal events in each fine bin is determined using a sideband-subtraction technique. 
The ratios of the acceptances and efficiencies, listed in Table [l| are different from unity because 
of small differences in the X(3872) and tp{2S) decay kinematics. 

Studies are performed to verify the description of the data by the simulations and to determine 
the systematic uncertainties. These are listed in Table |2]and described in the following. 

• Fit functions. The systematic uncertainty in the signal extraction from the invariant- 
mass spectrum is determined by variation of the fit parametrization independently 
for the X(3872) and tp{2S). Using a third-order Chebyshev polynomials for the back- 
grounds or the sum of a Gaussian and a Crystal Ball [24J function for the signal, 
variations of 1-2% are found. Fixing the X(3872) and tp{2S) mass difference to the 
PDG value il4| in the fit changes the result by less than 1%. 

• Muon-pair efficiency. Systematic uncertainties in muon efficiencies largely cancel in 
the cross section ratio measurement. Single-muon efficiencies are determined from 
]/ip events using a tag-and-probe technique on both the data and simulation |IT2| . 
The systematic uncertainty in the cross section ratio from this source is less than 1%. 

• Pion-pair efficiency. The systematic uncertainty in the efficiency for the reconstruction 
of the pion pair is determined by comparison of the measured and simulated event 
yields from ^(2S) — ?■ ]/xpn^n^ and ^{2S) — > ji^fi^ decays. After corrections for 
the branching fractions |1T4| and differences in the acceptance and efficiency for the 
muon pair, the ratio of event yields in the two decay channels differs from unity 




fbin 



rbin 



(3) 
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Table 2: Summary of the relative systematic uncertainties for Kfiduciai ^rid R. The variation over 
the Pj bins is given. The systematic sources common to both Rgduciai ^rid R are reported at the 
top, followed by those affecting only Rfidudal ^rid only R. 



Source 


Relative uncertainty (%) 


C^ommon to Rfj..^;^i and R 




Fit functions 


1-2 


Vr r / 


< 1 


e(7r+7r-) 


1-5 


Efficiency statistical precision 


1-3 


Specific to Raducml 




X(3872) pt spectrum 


2-5 


!/?(2S) Pj spectrum 


1-4 


Total systematic uncertainty in Kfiducial 


4-8 


Specific to R 




X(3872) pt spectrum 


1-11 


!/'(2S) Pt spectrum 


1-4 


m{n^n~) spectrum 


1-2 


Acceptance statistical precision 


1-3 


Total systematic uncertainty in R 


5-13 



because of different dipion reconstruction efficiencies. The more precisely measured 
value B{ip{2S) e+e^) |jl4l is used, instead of that for J/'(2S) — > fi^fi^, assuming 
lepton universality. Comparison of the simulation with the data reveals differences 
in dipion efficiency of 5% for pj < 15 GeV and 1% at higher transverse momentum. 

• Efficiency statistical precision. The efficiency uncertainties introduced by the statistical 
limitations of the simulated samples is less than 1% in general, rising to 3% for 30 < 
pi < 50 GeV. 

• X(3872) Pj spectrum. The dependence of the measurement on the transverse mo- 
mentum spectrum of the X(3872) is estimated by repeating the analysis with a sim- 
ulation including colour-octet contributions [21. 221- Simulations with and without 
colour-octet contributions lead to large variations of the pj spectra that are still com- 
patible with the data. The differences between these two cases, 2-5% on Rfiduciai ^rid 
5-6% on R, are taken as the systematic uncertainty. Variations of similar size are ob- 
tained when reweighting the simulated X(3872) pj spectrum to match the data. The 
uncertainty in the pj spectrum extracted from the data is also considered as a source 
of systematic uncertainty, which is added in quadrature. The uncertainty in Rfiduciai 
is found to be 2-5% . The rapidly changing acceptance as a function of transverse 
momentum makes the R measurement very sensitive to the pj spectrum, in particu- 
lar for low transverse momentum and for the px"iritegrated result. The uncertainty 
in R is 11% in the first pj bin and 1-7% elsewhere. 

• t/'(2S) Pj spectrum. For the xp{2S), the simulated pj spectrum is reweighted to match 
the distribution observed in data, and the efficiency and acceptance corrections are 
recalculated. The change in the cross section ratios, both for Rgdudai and R, is about 
4% in the lowest transverse momentum bin and 1-3% elsewhere. 

• m{n^TT^) spectrum. The dipion invariant-mass spectrum of the X (3872) ]/xpTT'^TT~ 
decay is extracted from the data, as described in Section |7| and compared with the 
expectations from the simulated samples. The dependence of the efficiency correc- 
tions on the dipion invariant mass is weak, and the systematic uncertainty in Rfiduciai 
is negligible. The dependence of R on the assumed invariant-mass spectrum of the 
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5 Measurement of the nonprompt fraction 



pion pair is estimated by reweighting the generated dipion invariant-mass spectrum 
to match the data. This leads to changes in the cross section ratio R of up to 2%. 

• Acceptance statistical precision. The uncertainty in the estimate of the dimuon and 
dipion acceptances owing to the statistical limitations of the simulated samples is 
1%, rising to 3% at high transverse momentum. 

The stability over time of the — ?► pi^j^^ yield relative to the ^{2S) — > ]/ipn^ W yield verifies 
that the muon and track selections used in the analysis are not sensitive to beam conditions or 
the amount of pileup. 

Adding all the systematic uncertainties in quadrature, a total systematic uncertainty of 4-8% 
in Rfiducial arid 5-13% in R is obtained. 

The cross section ratio is determined as a function of the transverse momentum of the ]/ipn^ n~ 
system. The results for both Rfiducial (Eq. and the fully acceptance-corrected R (Eq. ^) are 
listed in Table |3] and shown in Fig. |3] No significant dependence on transverse momentum is 
observed for either quantity. These results are obtained under the assumption that the X(3872) 
quantum numbers are J^'- = 1++, as favoured by existing data HSHHI, and no systematic 
uncertainty is assigned to cover other cases. 

In the simulations, unpolarized X(3872) and tp{2S) states are assumed. To evaluate the impact 
of other polarization scenarios, it is assumed that the X(3872) and the ]/xp from the X(3872) — > 
]/ipn^W decay have the same polarization. The polarization of the ]/'ip is varied in extreme 
scenarios, corresponding to fully longitudinal or fully transverse polarization in the helicity 
and Collins-Soper frames [25]. The same variations are performed separately for the ]/\p from 
X(3872) decays and for the tp{2S). The observed relative shifts of the cross section ratio R are 
listed in Tablejlj Small effects are found in scenarios where both the X(3872) and the ip (2S) have 
the same polarization. Assuming only one of the two states is unpolarized leads to variations 
of about 30% in the helicity frame and up to 20% in the Collins-Soper frame. Scenarios with 
transversely polarized X(3872) and longitudinally polarized tp{2S) give variations of up to 90% 
for the helicity frame and 30% for the Collins-Soper frame. In contrast, the fiducial cross section 
ratio Rfiducial is largely insensitive to polarization assumptions, showing maximal variations of 
4%. 

5 Measurement of the nonprompt fraction 

The relative contribution to the total X(3872) yield resulting from decays of B hadrons, often 
referred to as the nonprompt fraction, is determined from the decay lifetime distribution. The 
measurement is performed with the same ]/xp and tt^tt^ acceptance criteria presented above. 

Table 3: The ratios of the measured cross sections times branching fractions, Rgdudai arid R, as 
a function of the transverse momentum of the j/xpn^n^ system, together with their statistical 
and systematic uncertainties, respectively. For the first bin in transverse momentum, only the 
data from the period 2011a are included. 



Dataset 


Pt (GeV) 


Rfiducial 


R 


2011a 

2011a+b 

2011a+b 

2011a+b 

2011a+b 


10-13.5 
13.5-15 
15-18 
18-30 
30-50 


0.0767 ± 0.0082 ± 0.0059 
0.0728 ± 0.0076 ± 0.0044 
0.0724 ± 0.0056 ± 0.0042 
0.0611 ± 0.0042 ± 0.0025 
0.075 ± 0.012 ± 0.004 


0.0727 ± 0.0079 ± 0.0097 
0.0671 ± 0.0072 ± 0.0044 
0.0687 ± 0.0055 ± 0.0051 
0.0601 ± 0.0042 ± 0.0042 
0.078 ± 0.013 ± 0.004 


2011a+b 


10-50 


0.0694 ± 0.0029 ± 0.0036 


0.0656 ± 0.0029 ± 0.0065 
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Figure 3: Ratios of the X(3872) and ^{2S) cross sections times branching fractions, without 
(Rfiducial/ left) and with {R, right) acceptance corrections for the muon and pion pairs, as a 
function of pj. The inner error bars indicate the statistical uncertainty and the outer error bars 
represent the total uncertainty. The data points are placed at the centre of each pj bin. 



Table 4: Relative variations, in percent, of the integrated cross section ratio R for different 
X(3872) and ^(2S) polarization hypotheses: transversely (longitudinally) polarized ]/xp are 
denoted as CST (CSL) in the CoUins-Soper frame and HXT (HXL) in the helicity frame. Unpo- 
larized scenarios (labelled unpol) are also included. 



Polarization 


Relative 


Polarization 


Relative 


X(3872) 




shifts (%) 


X(3872) 


ip{2S) 


shifts (%) 


CST 


CSL 


-28 


CST 


unpol 


-8 


CSL 


CST 


+31 


CSL 


unpol 


+22 


HXT 


HXL 


+86 


HXT 


unpol 


+28 


HXL 


HXT 


-49 


HXL 


unpol 


-31 


CST 


CST 


-1 


unpol 


CST 


+8 


CSL 


CSL 


-5 


unpol 


CSL 


-22 


HXT 


HXT 


-6 


unpol 


HXT 


-27 


HXL 


HXL 


-1 


unpol 


HXL 


+25 



The "pseudo-proper" decay length ixy is defined in the plane transverse to the beam direc- 
tion as the distance between the vertex formed by the four tracks of the j/ipn^n^ system and 
the closest reconstructed primary vertex along the beam direction, corrected by the transverse 
Lorentz boost of the j/ipn^n^ candidate. An event sample enriched in X(3872) candidates 
from B decays is selected by requiring that ixy be larger than 100 jim. This selection retains 
about 80% of the nonprompt X(3872) candidates, while the contribution from prompt X(3872) 
is smaller than 0.1%, as determined from simulation. The simulated ixij distribution is veri- 
fied using the corresponding distribution from the t/'(2S) data sample. The nonprompt fraction 
is then obtained from the ratio between the signal yields in this B-hadron-enriched sample 
and the signal yields in the inclusive sample, after correction for the efficiencies of the decay- 
length-selection criteria, as determined from simulations of prompt and nonprompt X(3872) 
states. The signal yields are extracted from fits to the ]/xpn^n^ invariant-mass spectrum, as 
described in Section |4] In the fits to the B-hadron-enriched sample, the fit parameters for the 
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Figure 4: The J/t/^/r+Tr^ invariant-mass distribution in the X(3872) region for pj = 10-13.5 GeV 
(left) and 18-30 GeV (right) in a B-enriched sample. The lines represent the signal-plus- 
background fits (solid) and the background-only components (dashed). The ;^^/ndf of the 
fit is also reported. 



Table 5: Summary of systematic uncertainties in the X(3872) nonprompt fraction. 



Source 


Relative uncertainty (%) 


Vertex estimation 


1 


Background parametrization 


2-3 


Efficiency 


3-8 


Decay length resolution 


4 


Pileup 


2 


Total systematic uncertainty 


6-10 



mass and width are fixed to those determined from the full sample. Figure |4] shows examples 
of fitted invariant-mass distributions for the B-hadron-enriched sample. 

The measurement is found to be stable with respect to variations of the pseudo-proper-decay- 
length requirement between 50 and 250 jim. Application of the same method to ^{2S) candi- 
dates yields the same result as previously measured |T2| . In an alternative method, similar to 
the one used in Ref. [12], a two-dimensional fit to the invariant mass and the pseudo-proper 
decay length is performed. The £xy resolution is described by a function that depends on the 
uncertainty in the pseudo-proper-decay-length measurement, as determined event-by-event 
from the covariance matrices of the fits to the primary and secondary vertices. This function 
is obtained for signal and background, respectively, from the invariant-mass distribution after 
sideband subtraction [26| and from the sideband regions. These distributions are used to fix 
the lifetime parameters in the two-dimensional fit to correctly model the pseudo-proper-decay- 
length resolution. The validity of both methods has also been verified with simulated prompt 
and nonprompt X(3872) events for the signal, and combining from data the ]/xp and same-sign 
tracks into a j/xprcn candidate for the nonresonant background. Repeating the analyses on these 
samples, the nonprompt fractions are consistent with those used in the simulation. While both 
methods agree, the method utilizing the requirement on £xy is chosen since it has the smaller 
systematic uncertainty. 

Detailed studies are performed to determine the systematic uncertainties listed in Table |5] and 
described in the following. The selection of the primary vertex is modified by choosing the 
vertex with the smallest impact parameter along the beam direction for the X(3872) candidate. 
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Figure 5: Measured X(3872) nonprompt fraction, uncorrected for acceptance, as a function of 
Pj. The inner error bars indicate the statistical uncertainty and the outer error bars represent 
the total uncertainty. The data points are placed at the centre of each pj bin. 



Table 6: The X(3872) nonprompt fractions, not corrected for acceptance, as a function of the 
transverse momentum, together with their statistical and systematic uncertainties, respectively. 



Dataset 


Pj (GeV) 


X(3872) nonprompt fraction 


2011a 

2011a+b 

2011a+b 

2011a+b 

2011a+b 


10-13.5 
13.5-15 
15-18 
18-30 
30-50 


0.272 ± 0.057 ± 0.016 
0.182 ± 0.052 ± 0.013 
0.246 ± 0.043 ± 0.015 
0.297 ± 0.042 ± 0.021 
0.301 ± 0.097 ± 0.030 


2011a+b 


10-50 


0.263 ± 0.023 ± 0.016 



instead of the one closest to the four-track vertex along the beam direction. This variation 
changes the measured nonprompt fraction by 1%. The systematic uncertainties related to sig- 
nal extraction, determined by changing the background functions, are 2-3%. The difference 
between the reconstruction efficiency for prompt and nonprompt production, 8% for the high- 
est transverse momentum bin and 3-4% elsewhere, is taken as the systematic uncertainty. The 
uncertainty from the simulation of the pseudo-proper-decay-length resolution is estimated by 
comparing the ixy distribution from a simulated 1/^(25) sample with that from data. The change 
in the nonprompt fraction when relying on the £xy resolution from data is 4%. Finally, the sys- 
tematic uncertainty from the description of pileup events is evaluated from the dependence of 
the result on the number of primary vertices in the event and estimated to be 2%. From these 
estimates a total systematic uncertainty of 6-10% is obtained. 

The final results are listed in Table |6] and shown in Fig. [sjas a function of pj. The X(3872) 
nonprompt fraction reveals no significant dependence on transverse momentum and the inte- 
grated value is significantly smaller than that for the ^{2S) [12J. The results are obtained under 
the assumption that effects related to the X(3872) polarization cancel in the nonprompt fraction 
measurement, and therefore no systematic uncertainty is assigned for polarization effects. 
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6 Determination of the prompt X(3872) production cross section 

The cross section times branching fraction for prompt X(3872) production is determined from 
the measurement of the cross section ratio and the nonprompt fraction, described above, com- 
bined with a previous result of the prompt ^{2S) cross section fl2j. The latter measurement 
was performed using the tp{2S) — )■ decay mode and provides results as a function of 

transverse momentum up to 30GeV and for the rapidity range \y\ < 1.2. The prompt X(3872) 
cross section times branching fraction into ^/xpn^n^ is given by 



JX{3872) 



• R • (7' 



ip{2S) 



prompt 
V(2S) 



B{ip{2S) -^]/ip7T+n' 
~B{^p(2S)^]F^ 



where 'T'^ps)^' " '^(V'l^S) ) is the measured prompt J/'(2S) cross section times J/'(2S) 



branching fraction [12J, R is the cross section ratio reported in Sectionjij, and fx[2,9,7i) ^^'^ 
/i/i(2S) nonprompt fractions for X(3872) and 1/^(25), respectively. In the calculation, the 

branching fraction 5(i/?(2S) — > J/t/^/r+TT^) is taken from Ref. [14J, and 6(t/j(2S) — > is 
taken to be equal to the more precisely known B{^p{2S) — ?► e+e^) [14J. 

The corresponding differential cross section for prompt X(3872) production times the branch- 
ing fraction to J/ ^ tt^ as a function of transverse momentum, in the rapidity region < 1.2, 
is listed in Table [7| and shown in Fig.|6] No cancellation of systematic uncertainties is assumed 
in the combination. The main sources of systematic uncertainty are related to the measurement 
of the ratio R and the background lifetime fit in the measurement of the prompt 1/^(25) cross 
section l(T2ll . A calculation of the predicted differential cross section for prompt X(3872) produc- 
tion in pp collisions at = 7 TeV has been made using the NRQCD factorization formalism, 
assuming the X(3872) is formed from a cc pair with negligible relative momentum [llj. This 
calculation is normalized using Tevatron measurements [9,. 27} with the statistical uncertainty 
obtained from the experimental input data. The predictions from Ref. IITTI were modified by the 
authors to match the phase-space of the measurement presented in this paper. Comparisons of 
this prediction with the data, in Fig. |6| demonstrates that, while the shape is reasonably well 
described, the predicted cross section is much larger than observed in data. 

The integrated prompt X(3872) cross section times branching fraction for the kinematic region 
10 < Pt < 30GeV and \y\ < 1.2 is also determined. In this kinematic region, the ratio of 
cross section times branching fraction for X(3872) and ip{2S) is R = 0.0682 ± 0.0032 (stat.) ± 
0.0065 (syst.), and the nonprompt X(3872) fraction is 0.260 ± 0.024 (stat.) ± 0.016 (syst.). From 

Table 7: Prompt X(3872) differential cross section times branching fraction 5 (X (3872) — ?► 
]/\pn^ n^) as a function of transverse momentum of the ]/\pn^Ti^ system. The uncertainties 
shown are statistical and systematic, respectively. 



Pt (GeV) 


d£7-P™^PJ/dpT ■ fi(X(3872) ^ ]/ipn+TZ-) (nb/GeV) 


10-13.5 


0.211 ±0.034 ±0.035 


13.5-15 


0.081 ±0.013 ±0.010 


15-18 


0.0390 ± 0.0054 ± 0.0042 


15-18 


0.0390 ± 0.0054 ± 0.0042 


18-30 


0.0068 ± 0.0009 ± 0.0009 
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Figure 6: Measured differential cross section for prompt X(3872) production times branching 
fraction 5 (X (3872) — > ]/xp7i^7i^) as a function of px- The inner error bars indicate the statistical 
uncertainty and the outer error bars represent the total uncertainty. Predictions from a NRQCD 
model l|TT1l are shown by the solid line, with the dotted lines representing the uncertainty. The 
data points are placed where the value of the theoretical prediction is equal to its mean value 
over each bin, according to the prescription in |28j|. 

these results, the measured integrated cross section for prompt X(3872) production times branch- 
ing fraction is: 



£^prompt^pp ^ X (3872 )+ anything) ■ 6(X(3872) ^ ^/xpn+n-) = 1.06 ±0.11 (stat.) ±0.15 (syst.) nb. 

This result assumes that the X(3872) and 1/^(25) states are unpolarized. The NRQCD prediction 
for the prompt X(3872) cross section times branching fraction in the kinematic region of this 
analysis is 4.01 ± 0.88 nb |[Tll| , significantly above the measured value. 

7 Measurement of the 7z+7z~ invariant-mass distribution 

The decay properties of the X(3872) are further investigated with a measurement of the tt+tt^ 
invariant-mass distribution from X(3872) decays to ]/xpn^ . Here, the same event selection as 
described in Section|3]is applied. The event sample 2011a is used, with a transverse momentum 
threshold of 7 GeV for the muon pair, within the kinematic range 10 < px < 50 GeV and \y\ < 
1.25 for the J/j/^tt+tt^. In this sample, the X(3872) yield with the tt+tt^ invariant mass larger 
than 0.5 GeV is determined from a fit to the ]/tpn^n^ invariant-mass spectrum to be 6302 ± 
346, where the uncertainty is statistical only. The m(7r+7r^) > 0.5 GeV criterion is imposed to 
remove events with low efficiency owing to the requirement on the Q value of the decay. 

To extract the dipion invariant-mass spectrum from X(3872) decays, the event sample is di- 
vided into twelve intervals of dipion invariant mass in the range 0.5 < m{n^n^) < 0.78 GeV. 
In each interval, a maximum-likelihood fit to the ]/tpn^n^ invariant-mass distribution is per- 
formed, where the signal is modelled with a single Gaussian. The position and width of the 
X(3872) signal are fixed to the values obtained in the fit to the full sample, except for the last 
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m(J/\|/ [GeV] m(J/v ti-'ti") [GeV] 

Figure 7: Examples of the J/i/'tt+tt^ invariant-mass spectrum for the dipion invariant-mass 
intervals 0.67-0.69 GeV (left) and 0.73-0.75 GeV (right). The lines represent the signal-plus- 
background fit (solid) and the background-only component (dashed). 



interval, 0.765-0.78 GeV, where the mean and width of the Gaussian are left free to accommo- 
date possible distortions of the signal shape near the upper kinematic limit. The background 
shape in m{n^n^) intervals is different from the one for the entire m^n^n^) spectrum, and a 
third-order Chebyshev polynomial is used to model it, with the parameters left free in the fit. 
The ]/tp 71^71^ invariant-mass spectra for two of the n^n^ invariant-mass intervals are shown 
in Fig.[7| 

The X(3872) dipion invariant-mass distribution is extracted from the signal yields obtained 
from the fits to the data in each interval, after correction for detector acceptance and efficiencies, 
as estimated from the simulation. The resulting dipion invariant-mass spectrum, normalized 
to the total cross section in the interval 0.5 < m(7r+7r^) < 0.78 GeV, is presented in Fig. |8] The 
data are compared to X(3872) signal simulations with and without an intermediate in the 
]/ipTZ^n^ decay (generation details are described in Sectionjs]). The assumption of an interme- 
diate decay gives better agreement with the data, confirming previous measurements |l7ll20|. 

Detailed studies are performed to determine the systematic uncertainties. Scenarios with and 
without an intermediate p^ provide acceptance and efficiency corrections that are very sim- 
ilar. The impact on the acceptance correction from uncertainties in the X(3872) transverse- 
momentum spectrum is found by varying the simulated px spectra and generated pj distribu- 
tion to match the data. Variations of the corrected yields by 4-6% are observed and considered 
as a systematic uncertainty. 

The fits to the invariant-mass distributions are done with both free and fixed X(3872) mass 
and width. In addition, for modelling of the background in the higher dipion invariant-mass 
bins, a convolution of an exponential and an error functions is used, with a turn-on value con- 
strained to be close to the kinematic limit for each m ( 7r+ ) bin. These variations yield maxi- 
mal variations of the yields by 10-20%, and constitute the dominant systematic uncertainty in 
the measurement of the dipion invariant-mass distribution. 
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Figure 8: Dipion invariant-mass spectrum for X(3872) — ?► J/xpn^n^ decays corrected for ac- 
ceptance and efficiency. The distribution is normalized to unity by dividing by the total cross 
section for 0.5 < m ( 7r+ ) < 0.78 GeV. The inner error bars indicate the statistical uncertainty 
and the outer error bars represent the total uncertainty. The results are compared to results 
from EVTGEN with (solid line) and without (dotted line) an intermediate decay. 

8 Summary 

The X(3872) production cross section has been measured in pp collision at a/s = 7TeV, with 
data collected by the CMS experiment, corresponding to an integrated luminosity of 4.8 fb^^. 
The measurement makes use of the decays of the X(3872) and ^{2S) states into ]/tpK^n^ , with 
subsequent decay of the ]/xp into two muons. The ratio of the inclusive cross section times 
branching fraction of the X(3872) and tp[2S) in the kinematic region 10 < px < 50 GeV and 
lyl < 1.2 is R = 0.0656 ± 0.0029 (stat.) ± 0.0065 (syst.). When restricted to the measured phase- 
space of the muon and pion pairs, the ratio is Rfidudai = 0.0694 ± 0.0029 (stat.) ± 0.0036 (syst.). 
These ratios show no significant dependence on the transverse momentum of the ^/xpn^n^ 
system. The results have been obtained with the assumption that the X(3872) has quantum 
numbers }^'~ = 1++ and that both the X(3872) and the ip{2S) are unpolarized. Variations of the 
results for different polarization assumptions have also been reported. The fraction of X(3872) 
originating from B-hadron decays is 0.263 ± 0.023 (stat.) ± 0.016 (syst.), again assuming the 
X(3872) is unpolarized. No significant dependence on transverse momentum is found. From 
these measurements, the cross section for prompt X(3872) production times branching frac- 
tion into ^/tpn^n^ has been extracted, using a previous CMS measurement of the cross section 
for prompt xp{2S) production. A value of aP'^^'P^pp X(3872) + anything) ■ 6(X(3872) 
]/tpn^7i^) = 1.06 ± 0.11 (stat.) ± 0.15 (syst.) nb is found for the kinematic range 10 < pj < 
30 GeV and \y\ < 1.2. This result is also made under the assumption that the X(3872) and 
i/?(2S) states are unpolarized. The NRQCD predictions for prompt X(3872) production at the 
LHC significantly exceed the measured value, while the pj dependence is reasonably well de- 
scribed. The measured dipion mass spectrum for X(3872) — ?► ]/xpn^n^ clearly favours the 
presence of an intermediate state. 
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